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We suggest an effective field theory for disorderd conductors, which describes quantum kinetics 
of ballistically propagating electrons. This theory contains non-linear cr-model (l[ as its long wave 
limit. 

- - -, 1. The non-linear cr-model is proven to be a useful tool in the description of various properties 
Q\ ", of disordered conductors. Any property, such as conductivity, averaged over different realizations of 



the random potential can be presented in this model as a statistical average with the free energy 
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3 ■ F = — drstr[D(VQY + 2iujAQ), Z = / VQe~ b (1) 

^—5 ' 8 J J 

Q 2 =l 

The functional integral is taken over the 8x8 super-matrix Q(r) which is subjected to the constraint 
Q 2 — 1. Here and below we use the super-matrix version [|TJ of the nonlinear a model. 
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■ This discription is valid under the following two conditions: 

o\ : 

i) The Fermi wave length = fi/pp is much smaller than the mean free path /, i.e. ppl/fi ^> 1. 

■ ii) The typical wave vector q of the super-matrix fluctuations is smaller than 1/Z, i.e. ql -C 1. 

These conditions mean that (i) the semi-classical description is applicable to the electrons with 
the Fermi energy, and (ii) their motion is described by the diffusion equation. There are physical 
situations when the condition (i) is fulfilled, while the condition (ii) is not and electrons propagate 
ballistically. This happens, for example in a metallic grain with a diffusive boundary scattering if 
q ! the bulk mean free path I is much larger than the grain size L, i.e. / ^> L. 

j> ! In this letter we present a generalized version of the model (P whose validity is no longer restricted 
^ by condition (ii). The generalized partition function correctly accounts for the fluctuations with wave 
vectors q ~ l/l and therefore can be used for the description of systems with ballistic electron motion. 

_ _ i 

We begin with a general expression for the free energy which is obtained after averaging over the 
random potential, the Hubbard- Stratonovich decomposition of the quartic form and integration over 
the electron degrees of freedom (see QXJ] for details and notations). 

F = ~stv\n[-iK] + ^y s trQ 2 (r)dr, Z = J VQe~ F , (2a) 

A = £-#o + -A + -Q, H = HW)2 (2b) 
2 2r 2m V ; 

This expression appears at a preliminary stage in the derivation of Eq. (P and the supermatrix Q 
is not yet restricted by the constraint Q 2 = 1. 

Equation (|2a|), in principle, could have served as a required generalisation of the free energy (jl]). 
However, it is too detailed being valid for the super-matrices Q fluctuating with arbitrary wave 
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vectors q. It will be simplified in order to describe the small q fluctuations only (q <C pp/h). The 
first step in the simplification is the same as in the derivation of the quantum kinetic equation in 
the Keldysh approach (see, for example, M). 
2. The Green function G(r, r'\Q) of the operator K obeys the equations 



E - floOO + f A + J-Q(r) 



E-H (r') G(r,r'\Q)+G(v,v'\Q) 



G(r,r'|Q) = i5(r - r') 
= z5(r — r') 



-A H Q(r' 
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Subtracting Eq. ([T5|) from Eq. fl3"a| ) and going to the Wigner representation 



G(r,r') = (dp)G( 



~ ,r + r' 



p) e ip(r - r,) 

we can find after the integration over the modulus of the momentum p an equation for 
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g n (r) = - d£G(r,n±- 



This equation can be presented in the form 

%i(r) 



2t>^n- 
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(3a) 
(3b) 

(4) 
(5) 
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which resembles the quantum kinetic equation in the Eilenberger form ||. The matrix g n (r) in this 
equation has the meaning of distribution function at a coordinate r and momentum p = n • pp. 

Being linear, Eq. (|| does not define g n uniquely and must be supplied with the normalisation 
condition [01 



#n = 1 ; tr# n = 0. 

The matrix Q(r) is invariant with respect to the charge conjugation 

Q = CQ T C T = Q, 

where C is a certain matrix (see Q), C T C = 1. Taking the charge conjugate of Eq 
Eq (|), we see that G(r, r') obeys Eq (|3t|). Therefore 

G(r,r') = G(r',r), d(r, p) = G(r, -p), ^ n (r) = p_ n (r). 



(7) 



and using 



(9) 



Thus, Eq. (P) with the normalisation condition (0) and the symmetries (||]) is a long wave limit of 
Eqs. (0). Our goal is to perform analogous simplification of the free energy (|2a|). 

3. An intermediate step is finding a functional $, which reaches its extrema for solutions of Eq. @. 
This equation resembles the equation of motion of a magnetic moment M in external magnetic field 



B : 



dM 



[M x B] 
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The action for this problem has the form (see, for instance, 0) 
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A= I dt'BM(t') + J dt' J duM- 



dM dM 

dt du 



where the function M(t, u) is introduced as 

M(t,0) = M ; M(t,l) = M(t). 



(11) 
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The second term in Eq ( |TTD does not depend upon the choice of M and values of M(t, u) for 
< u < 1, provided M(0) = M(t). 
Following this analogy we present $ in the form 
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$ = J drstr 
(9(r)) - 



-Q{v)-iuA){g{v)) 



dQ r 



■in 



dQ n r 1 



W{# n } = J dr J J du sir g n (r, 
# n (r, 0)=A; g a (r, 1) = g a (r). 



u 



dg n dg Y 
n- 



du dr 



(13a) 

(13b) 

(13c) 
(13d) 



The functional derivative 5Q/5g n must be taken with constraint (|7|) which guaranties that g n Sg a + 
figngn = and an arbitrary variation 5g n has the form Sg n = [g n , a n }. As a result 



5$ = dr 



dQ n 

47T 



str 



-Q(r) - iuA,g n 

T 



where 



SW = 4 / dr 



dQ n 



, / dg n 
str n— — a 



or J 



(14a) 



(14b) 



Thus, Eq. ( |13a| ) gives the required functional. 

4. Now we are prepared to show that in the limit / ^> \p the partition function (pa]) reduces to 
the form 



Vg n {r)e- 



F= — 
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J rfrstr [iuA(g(r)) - ^( 9 (r)) 2 } - y ■ MU 
W{9n} = J dr J du str g a (r, ■ 



dg n n 9g u 
du ' <9r 



Indeed, the following indentity is valid 



Zi{Q} = exp 



- str ln(— iK) 



Vg n (r) exp 



7TZ/ 



^ 2 {g}. 



(15a) 

(15b) 
(15c) 

(16) 



The free energy 7rz/$/4 in the partition function Z 2 has a deep minimum for g n (r) equal to g^\r\Q) 
which is the solution of Eq. (0). With the saddle point precision 
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SZ 2 {Q} _ 



^ J 0/(r)) exp 
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(17) 



On the other hand 



SZ^Q} Z^Q} 7tz/ (0) 

= ~^ G(r ' r) = ^ {9n } ' Zi{Q} - (18) 

Thus, the functionals Z 12 {Q} obey indentical equations. Since Zi{A} = Z 2 {A} = 1 the indentity 
(16) is proven. 

Substituting Eq. ([H]) into Eq ( pa] ) and taking the Gaussian integral over Q, we arrive at the final 
expression (|i~5|). 

5. For small gradients, the free energy (pj|) reduces to the standard cr-model (|l]). To show this we 
expand the matrix g n into the sum over sperical functions Yl,m(h) 

oo L 

9n(r) = J2J2 9l,m( t ) ■ Y L,Al( n ) 

L=0 M=-L 



and note that only zero and first harmonics contribute to the functional integral (|i~5|) : 

QJ 2 



9n = Q(r) + J(r) ■ n 
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The constraint g 2 = 1 now reads 



qj + jg = o. 



(19) 



(20) 



Substituting the Eq. fll9|) into Eqs. ([15|) and using conditions (^0]) we obtain the partition function 
in the form 

J 2 



Z = JvqJ VJe- F(Q ' J \ F(Q,J) = ™J drstv{iuAQ+ — - ^(VQ)QJ} 



(21) 



After the Gaussian integration over J in Eq (M) we arrive, finally, at Eq (|T|). 

6. Equations ( |15|) can be generalised in order to describe the ballistic motion in the presence of 
external fields. In a general case the electron is described by the classical hamiltonian H(pi,Xi) and 
the kinetic equation ([]) has the form (see 0): 

/iuA Q 
A~2 2t, 



{#(x,p),#(a;,p)} 
where {if, g} denotes the Poisson brackets 



,g(p,x) 



(22) 



{H{x,p),g(x,p)} 



dH dg OH dg 



dpi dxi dxi dpi 

Equation (|2"2] ) is still the first order differential equation and the generalisation of expression (|15b| ) 
for the free energy has the form 



F = — [ dxidpibiE — H(p, x)) str \ iuAg — — - [ dug(x,p, 
4 7 2r 2 Jo 



u 



;>•» 



(23) 
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where 



dp>5(E-H(p>,x))g(x } p'). 



7. As an application of Eq. (p3|) , let us consider the derivation of the Pruisken action || for a 
two-dimensional electron gas in a perpendicular magnetic field B. To simplify the treatment, we 
consider only the case of classically weak field 

eB 

n c r < 1; tt c = , (24) 

mc 

when there is no Landau quantisation and the density of states v is a constant. Nevertheless, we 
take into account that in the presence of magnetic field the symmetry of g-matrix is reduced, and g 
belongs to the unitary ensemble. The Poisson brackets in magnetic field are 



{H,g} = v F n 



dg n 
dr 



a 



n x 



and the free energy (|23| ) has the following form 

■nv f [ (g) 2 If 1 I 

F=— J drstr <iuA(g) - — - -J du(g(x,p,u) 



dgn 
dn 



dg dg 

— ,v F n— + ll c 

on dr 



n x 



dg_ 
dn 



(25) 



(26) 



In the diffusive limit the expansion ([19]) can be used, which leads to the following expression for the 
free energy as a functional of Q and J; 

— T2 O 

Z = jvQ j PJe- F(Q ' J) , F(Q, J) = ^ J dr stv{iuAQ + — - y (VQ)QJ -yQ[JxJ]} (27) 

The last term in the free energy fl2"7] ) does not vanish because the components of the matrix J do 
not commute. Under the conditions ([24]), the Gaussian integration over J may be performed, with 
the vector product in Eq. (p7] ) as a perturbation, to yeild, finally, the free energy in the form 



8e 2 



drstr (a xx {VQ) 2 + 2a xy Q[V x Q,V y Q}) 



where 



On 



e 2 vD, 



a. 



(28) 



(29) 



8. There is a topological question, related to the W-term in the free energy fll5|): is it always 
possible to construct the functional W{g}, whose variation is given by Eq. (|14b|) ? The prescription 
( |18c| ) gives the W-term for the functions g{r), which are close to (?o( r ) = A. The question is whether 
such a functional can be defined globaly. 

The answer depends upon the topology of the constant energy surface H (r, p) = E in the phase 
space {xi,pi}. For the cases of billiards and space dimension d > 1 the functional W does exist. 

For a one-dimensional system W can only be found as a multivalued functional, just as the action 
(pT|). This causes no trouble, provided nhuvp is an integer. This integer exactly equal to the 
wave-guide channel number in the wire. 
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An accurate mathematical formulation and the proof of these statesments will be presented else- 
where. 

9. So far, we have considered only the systems with finite amount of disorder. One can see, however, 
that the expression (|23"D remains meaningful even as r — > oo. Therefore, we expect that the free 
energy = F(t — > oo) describes a clean system with the Hamiltonian H. As a consequence, the 
partition function = j Vg exp(— Foo) with the proper source terms gives the level statistics. 

In the low-frequency limit {u — > 0) only the zero-mode g°{r,p) such that {H, g } = contributes 
to Z^. There are two possibilities: 

i) The hamiltonian system under cosideration is integrable and there exists a set of integrals of 
motion {Ii, . . . I n }, {H,I k } = 0. Under this condition the energy levels are characterised 
by the eigenvalues of {Ii, . . . I n } and do not repel each other. Therefore the level statistics is 
Poissonian. 

ii) The classical dynamics is chaotic and the only integral of motion is energy. In this case the 
zero-mode is constant in the phase space and Z^ is reduced to the form 

Zoo= J £><?exp {-^ str(A<7)) (30) 

which leads to the Wigner-Dyson (WD) level statistics []]]. 

In the chaotic case deviations form the WD statistics occur for the frequencies larger than the 
inverse time of flight through the system. These deviations are described by the small fluctations 
of g about several stationary points Aj, similar to what has been recently shown by Andreev and 
Altshuler (AA) for diffusive systems 0. In complete agreement with a general AA-conjecture, the 
deviation from the WD statistics is described by the determinat of some operator. It follows from 
our consideration that this is the Liovillean operator 

- _ dH d dH d 
dp dx dx dp 

10. In conclusion, we would like to emphasize that the theory presented here contains the diffusive 
cr-model as a limiting case and supplies it with the physically motivated regularisation of the infinities 
at short distances. 

We are greatful to B.D.Simons for his question, which inspired us to derive Eq (|28|). 
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